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Abstract 

We derive numerical tight-binding hopping parameters to describe 
conduction-band dispersion in arbitrary orientational phases of A3 Ceo- The 
parameters are obtained by direct Fourier inversion of the spectra from self- 
consistent electronic-structure calculations for K3C60 using the local-density 
approximation, including the effects of orientational dependence. Using the 
new parameters, we revisit several earlier investigations of the orientational 
ordering in A^Cqq; some of the earlier results are substantiated, while others 
are slightly modified. 
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The electronic structure of the fullerene and fulleride sohds has been treated at vari- 
ous levels of theoretical sophistication, including simple Hiickel calculations for the isolated 
molecule, fully self-consistent density-functional methods using the local-density approxi- 
mation (LDA) for the intercalated solids, and most recently, quasiparticle-corrected GW 
methods for solid Cgo [JT]^]- For the A^Cqq family (x=0,3,4,6, A=heavy alkali), a single 
common picture emerges: Ceo molecular states are broadened into band states in the solid, 
with bandwidths on the order of ~1 eV or less, and the alkali s-electrons populate the 
Ceo-derived band states. The small dispersion, together with the one-to-one correspondence 
between molecular levels and bands, implies a nearly negligible redistribution of electronic 
charge density on forming the condensed phase (the alkah ionization aside). This suggests 
that a tight-binding (TB) Hamiltonian might provide an accurate description of dispersion 
in this family, that TB parameters should be transferable among similar structures, and 
that interactions may be limited to nearest neighbors. 

We focus on the superconducting A^Cqq phase because the nature of fullerene molecular 
orientational order in this phase is not well understood, having been variously described as 
merohedrally disordered with Fm3m symmetry |Q; orient at ionally disordered with lower 
symmetry |]^; and partially ordered on short length scales Important theoretical work 
by Gelfand and Lu and by Satpathy et al. |^ demonstrated that complete merohedral 
disorder has a profound effect on the conduction-band spectral density and on ac conduc- 
tivity, effectively washing out the sharp structure that normally results from a well-defined 
band structure. Deshpande et al. subsequently developed an effective-medium theory 
to describe the effects of disorder in a Green's function formalism, essentially confirming 



the results obtained numerically by Gelfand and Lu. Recently, Mele and Erwin [|T0| used 
both formal and numerical methods to show that although disorder strongly scatters electron 
states and substantially redistributes spectral density, states near the Fermi energy are char- 
acterized by an effective Bloch wave vector and a mean free path of order 2-3 near-neighbor 
spacings. Consequently, a relatively well-defined "Fermi surface" may be considered for 
these states even in the presence of random (merohedral) disorder. In all of these studies, 



the Hamiltonian of Ref. was used, for which hopping parameters were calculated from Ceo 
interactions at the level of Hiickel theory. In this Report, we provide an explicit parametrized 
Hamiltonian for A^Cqq, with TB parameters derived directly from first-principles LDA cal- 
culations of dispersion in A^Cqq, with the effects of orientational disorder fully accounted 
for. 

The structure and symmetry of the nearest-neighbor hopping Hamiltonian have been 
discussed in detail by Yildirim et al. ||Tl|] , and we will adopt their notation for the discussion 
here. The Hamiltonian describing conduction-band dispersion is 



R,r Q/3 



(1) 



where a,j3 = x,y,z is the orbital polarization, ctr = ±1 is an Ising variable giving the 
orientation of the molecule at R, and r = (a/2)(110), etc., gives the direction of the hop. 
The hopping matrix has four possible values, and can be conveniently decomposed as: 



(2) 



For the hopping direction tq = (a/2)(110), Yildirim et al. showed from symmetry consider- 
ations that the t^^^ take the form 



ABO 
B A 
C 



X Y 
-Y -X 




X -F 
F -X 




DEO 
EDO 
OOF 



(3) 



A different, equivalent linear combination of these eight constants was used by Gelfand and 
Lu, who defined hopping matrices between "parallel" (like orientation) and "perpendicular" 
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(unlike orientation) molecules, within a basis that rotates with the molecules. These two 
matrices have the form 



(4) 



To fit the constants in Eqs. and (Q) to self-consistent LDA spectra, we first consider 
a supercell of the A^Cqo fee Bravais lattice, with a basis of Ceo molecules with orientations 
c"Ri {i = If the orientations are not all identical, then all four terms from Eq. (0) 

will contribute to H. By calculating the LDA eigenvalues, E^^^, throughout the zone, 
optimal values (in the least-squares sense) for all of the constants in Tables I and II can be 
determined. We start by assuming we know the eigenvalues, -Ejif , and eigenvectors, |nk), 
of Ti. in the TB basis for a particular set of hopping constants. Then 7i(k) can be written 



as 



^(k)=E^SfW(nk|. 



(5) 



From Eq. (|T]), the matrix elements of ?i(k) in the TB basis are also given explicitly in terms 
of the hopping parameters: 



[R, a\H{k)\R + T, P) = Eta/3(r; aR, aR+r) exp(ik ■ 



(6) 



Eqs. (I) and (|^) are equivalent formulations. We exploit this by first constructing 7Y(k) 
according to Eq. (^ but with the substitution E'^^ E^^^. The resulting Hamiltonian 
matrix is Fourier inverted according to Eq. (^ to yield optimal hopping parameters in the 
least-squares sense. In practice, the solution is found iteratively, with the hopping constants 
from Ref. |^ providing the seed eigenvectors. 

The supercell we choose is the smallest possible, containing two Cgo molecules (with 
different orientations) and six K ions. The resulting crystal structure has D^h symmetry, 
and can be described as sheets of Ceo molecules alternating along the (001) direction between 



the "A" and "5" orientations; this corresponds to the ^2-62 structure of Yildirim et al. [|1T 
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The computational methods for solving the Kohn-Sham equations are identical to those used 



in previous studies of KgCgo, K3C60, and K4C60 fT^l- '^^^ charge density and potential are 
completely general and without shape approximation, and were iterated to self-consistency 
using the F point. Bloch basis functions are linear combinations of occupied and unoccupied 
atomic orbitals for potassium and carbon; these are in turn expanded on a set of gaussian 
functions. This provides a compact basis, and allows us to perform accurate all-electron 
calculations — without the need for pseudopotentials — in which core, valence, and conduction 
states are treated on equal footing. The Ceperly-Alder exchange-correlation functional was 
used. For the inversion of Eq. (^), we used a mesh of 40 equally spaced k-points in the 
irreducible wedge of the folded zone. 

The hopping constants obtained by fitting to the LDA spectrum (denoted TB/LDA) 
are listed in the first row of Tables I and II; the constants obtained by Gelfand and Lu 
(denoted TB/GL) are listed in the second row of Tables I and II. The TB/LDA values 
follow the qualitative trends in the TB/GL values, but substantial differences are apparent. 
To clarify the significance of these differences, we have used the Hamiltonian of Eq. (|T|) to 
compute spectra for three different structures: the orientationally ordered crystal (Fig. 1); 
the A2B2 supercell used to extract the hopping parameters (Fig. 2); and an ensemble of 
27-molecule supercells with quenched disorder, which simulates the merohedrally disordered 
structure originally proposed by Stephens (Fig. 3). For the first two structures, we also 
compare the TB spectra to the full self-consistent LDA spectra. We only show densities-of- 
states (DOS), since differences are still evident in these k-integrated spectra. 

Not surprisingly, for the orientationally ordered crystal, the TB/LDA parameter set 
reproduces the full LDA calculation more faithfully than the TB/GL set. More interesting 
is the DOS (Fig. 1, middle panel, light curve) using the hopping constants fit to the LDA 
spectrum from the single-orientation fee crystal. The close resemblance to the spectrum 
obtained by fitting to the LDA spectrum of the A2B2 structure illustrates the transferability 
of a particular hopping Hamiltonian from one crystal structure to another. This is the 
numerical justification for using hopping constants obtained from the A2B2 structure to 



compute spectra for supercells simulating merohedral disorder. Interestingly, the ensemble 
average of these supercell spectra (Fig. 3) shows only small differences arising from the two 
parameter sets. The small dip at 0.22 eV in the TB/GL data is absent in the TB/LDA 
spectrum, and the band edges of the latter show more gradual tailing than the former; 
otherwise the spectra are very similar. 

Deshpande et al. @ have discussed the fact that the magnitudes and signs of the constants 
in Tables I and II depend on the choice of local basis sets for two neighboring molecules. They 
defined a simple gauge-invariant quantity by considering the spectrum of the Hamiltonian for 
hopping around the smallest possible closed path, i.e. forward and backward hopping across 
a single "bond". There are two possibilities — hopping between like {AA) and unlike {AB) 
orientations — which the authors showed lead to numerically similar spectra, as shown in the 
second row of Table III. If it were true that hopping proceeded via nearest-neighbor atomic 
sites only, then the structure of the Ceo molecule and the crystallographic alignment of the 
two Cgo orientations would lead to exactly degenerate single-bond spectra. The deviation of 
the values in Table III from this equality is a direct measure of the importance of hopping 
between second-neighbor (and higher) atomic sites. In the first row of Table III, we have 
recomputed the single-bond spectra based on the new hopping constants. Again, changes 
are evident (especially for the smallest eigenvalue), but the approximate degeneracy is still 
observed. In particular, the largest AA and AB eigenvalues — which dominate the spectral 
features for all such retraceable paths — are nearly equal for both sets of constants. This 
justifies numerically the treatment of Ref. [^], in which all sites belonging to retraceable 
paths in the effective medium were assigned the same orientation. 

Finally, we reconsider the results of Yildirim et al, who have made a systematic analy- 
sis of the ground-state energetics of various orientational phases of A3C60 [|n|]- The authors 
considered the Ising-like nature of the orientational ordering problem, and showed that while 
the direct Coulomb interaction between molecules favors a "ferromagnetic" ordering, the in- 
direct kinetic energy contribution from the conduction electrons favors an antiferromagnetic 
alignment of nearest neighbors. The orientational dependence of the Coulomb interaction 
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was analyzed in detail in Ref. [0; this analysis was based directly on a multipole expansion 
of the electronic charge density calculated within LDA, and so remains unaffected by the 
present work. The kinetic energy contribution was analyzed in Ref. |]1T[ by calculating the 
energy of the filled Fermi sea, Eke = 2 J^nk ^nkfnk, for three different orientational phases: 
the orientationally ordered fee crystal, denoted by A4; the A2B2 structure discussed above; 
and a cubic (Pm3) structure denoted ^3!?. In Table IV, we have recomputed Eke for 
these three structures, using our TB/LDA hopping constants, as well the TB/GL constants 



originally used in Ref. |]lT 



Several aspects of the results in Table IV deserve comment: (1) Both sets of hopping 
parameters favor the non- ferromagnetic A^B and A2B2 configurations; (2) the TB/LDA set 
predicts the non-ferromagnetic configurations to be considerably more favored than does 
the TB/GL set; (3) the TB/LDA set predicts the energy ordering of the A^B and A2B2 
structures to be opposite to the prediction of the TB/GL set; (4) the ratio of relative 
energies, A{A2B2) / AlA^B) , has the value 1.06 for the TB/LDA set, considerably closer to 
the ideal nearest-neighbor Ising value of 4/3 than is the TB/GL value of 0.81. This last point 
suggests that the TB/LDA kinetic energy contribution to the orientational potential is also 
well modeled by a nearest-neighbor Ising model. When both Coulomb and kinetic energies 
are considered, the new parameter set implies that: (5) Non-ferromagnetic configurations 
are still favored so long as the kinetic energy scale, t, and the Ising interaction strength, 
Jd, satisfy the inequality t > —1.86 Jo- This condition is easily satisfied for all values of Jd 
considered in Ref. |jll| (note that the ferromagnetic Coulomb interaction implies Jn < 0). 
(6) Finally, the A^B and A2B2 structures may in principle be degenerate if t=— S.OOJ^^; the 
LDA value for t is smaller than this, but not unreasonably so. In summary, the new TB/LDA 
results suggest that the total orientational interactions is in fact reasonably well described by 
a nearest-neighbor antiferromagnetic Ising model, provided that current estimates of direct 
and indirect energy scales are meaningful. 
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TABLES 



TABLE L Values for the hopping constants appearing in Eq. (^). To allow comparison between 
the two sets of parameters, we set t=0.0146 eV, which equates the bandwidth of the Gelfand/Lu 
spectrum with the LDA bandwidth. 
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LDA fit 
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TABLE III. The three positive eigenvalues of the Hamiltonian for hopping back and forth 
across a single bond, for like (AA) and unlike {AB) orientations. 
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TABLE IV. Kinetic energy of the conduction electron states in three different orientational 
phases, in units of t. In the bottom panel the energies are given relative to the ordered A4 phase. 
Note the overall increase in magnitude of both relative energies in the LDA results, as well as the 
reversal of the ordering between the LDA and Gelfand/Lu results. 



St.mrt.iiT'P 




TCinptiP PTif*rp"v /molppnlp ftl 




LDA fit 


Gelfand/Lu 




-25.12 


-24.13 


A3B 


-28.18 


-28.08 


A2B2 


-28.43 


-27.53 


A3B - ^4 


-4.05 


-2.94 


^2^52 - ^4 


-4.30 


-2.39 



11 



FIGURES 

FIG. 1. Spectra for orientationally ordered fee A^Cqq, as calculated by self-consistent lo- 
cal-density methods (top), the tight-binding parameters fit to LDA spectra (middle), and the 
parameter set introduced by Gelfand and Lu (bottom). In the middle panel, spectra are shown 
for TB parameters fitted to LDA spectra of two different orientational phases, illustrating the 
transferability of these parameters. The dotted line is the Fermi energy. 

FIG. 2. Spectra for orientationally modulated A2B2 structure described in the text. Panels are 
labeled as in Fig. 1. 

FIG. 3. Spectra for an ensemble average of 27-molecule supercells with quenched disorder, 
using two different hopping-parameter sets. 
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